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ABSTRACT 


The damping of the ion-acoustic waves by Coulomb collisions 
is studied by using the Fokker-Planck equation of Rosenbluth et al. 
for both the species constituting the plasma, namely the electrons 
and the ions. In a plasma with weak collisions and in the absence 
of any external field, one finds that irrespective of the ratio of the 
ion temperature to the electron temperature, T, the characteristic 
frequency gets affected only by the electron-ion collisions. How- 
ever, as far as the collisional damping is concerned, both the 
ion-ion and the electron-ion collisions play a somewhat equally im- 
portant role; the electron-electron and the ion-electron contribu- 
tions are negligible compared to the other two. The damping in- 
creases with increase in the collision frequency but decreases with 
increase in T. 
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I. INTRODUCTION 


Ion acoustic waves, which were first predicted by Tonks and Langmuir 1 
using the fluid analysis, have been studied by a number of authors on the 
basis of collisionless Boltzmann equation. The collisionless theory shows that 
when the electrons and the ions have the same temperatures, these waves are 
very heavily Landau-damped. Experimentally, the ion waves were first observed 
by Revans® and their Landau-damping was measured by Wong et al.^ Wong 
et al. studied the space damping rather than the time damping in Cesium and 
Potassium and observed that the damping constant depends on the magnitude 
of the ion drift. 

The collisional damping of the ion waves was studied by Bhadra and Varma ; 8 
they used Krook’s model 9 to describe the ion-ion collisions. Kulsrud and 
Shen 10 used a slightly more realistic model for the collisional process; in their 
model, the ions were described by the Fokker-Planck equation of Rosenbluth 
et al ., 11 however, the electrons were treated by a fluid equation with the further 
assumption that the electrons were isothermal. This assumption of electrons 
being isothermal is equivalent to considering the Vlasov equation for the 
electrons, thereby neglecting the electron-electron and the electron-ion col- 
lisions. We remove this restriction on the electrons and treat both the 
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electrons and the ions by the corresponding Fokker-Planck equations. We in- 
deed find that the electron-ion contribution to the collision term can not be 
neglected as compared to the ion-ion contribution. 

Similar calculations for high-frequency plasma waves have been done by 
Comisarl2 and Buti. 1 **-!^ The effect of strong collisions on the low-frequency 
electrostatic plasma oscillations with and without an external magnetic field 
has been studied by Kuckes^ by using fluid equations. Following Comisar 
and Buti, we solve the Fokker-Planck equation on the assumption that the 
Coulomb collisions are weak i.e., car » 1, where oo is the characteristic fre- 
quency of the wave and r is the mean collision time for the collision process 
under consideration. We find that in the absence of any external field, irrespec- 
tive of the ratio of the ion temperature to the electron temperature, the char- 
acteristic frequency gets affected only by the electron-ion collisions; however, 
as far as the collisional damping is concerned, both the electron-ion and the ion- 
ion collisions play a role: the electron-electron and the ion-electron contributions 
are negligible compared to the ether two. 

II. GENERAL THEORY 

Consider an unbounded fully ionized hot plasma consisting only of electron 
and ions without any external field. In equilibrium both the electrons and the 
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ions obey Maxwellian distribution of velocities i.e., 


f o a (v) - ( 277V a 2 ) 3/2 e 


-v 2 /(2V a 2 ) 


(1) 


where the subscript a stands either for the electron or for the ion and V a 2 = KT a /m a 
with m e = m and m i = M. Nf oa is the equilibrium distribution function for 
species a. For small perturbations, the linearized Fokker-Planck equation 
for longitudinal oscillations is given by 


~dt + v 


at "f. Ho, 

dx + m a d\ 



( 2 ) 


where f a (x, v, t)is the perturbed distribution function and (df a/ /<9t) c takes 
into account the collisions between like and unlike particles and is represented 
by 



_d_ 

dv 



d 2 _ 
dv dv 


Tr 


/ A A \ 

\““/l 


-i- c 


/a a\ 1 

\ / O aj 


\ - / 


with 
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( 6 ) 


(aa) 3 - Nr a 2^ 5757 f dv ' f j (v '> 1 Y ~ V 'I 

J=e,i J 

and 

( AA )oa = NF a Y'. f dV ' f oJ O') U-V'l • 

J=e, i 

In these equations T a = {^n e a 2 /m a 2 ) ■in (477 NX. D 3 ) with ^ D 2 = KT av/ /(47rNe 2 ) ; T av is 
the average of the electron and the ion temperatures. The perturbed electric 
field in Eq. (2) is given by the Poisson’s equation 


div E 



( 8 ) 


To solve the pair of coupled equations (2) and (8), we take the Fourier trans- 
forms in space and Laplace transforms in time of all the perturbed quantities; 
Eq. (2) then becomes 


Ne a af oa ( df »\ 

(s + i k • v) f a ' (k,v, s) - g a (k, v) + — E k ' • = ^7t/ c ’ ( 9 ) 

where f a and E k ' are the Fourier- Laplace transforms of f a and E and g a (k, v) 
is the Fourier transform of the initial perturbation in the distribution function. 
To ensure the convergence of the integrals, we can have s in Eq. (9) such that 
Re s > 0 . On further taking the Fourier transforms in velocity space, thus 
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defining 


F (k, a, s) 


J"dv 


- ia • v 


f ' (k, v, s) 


(10) 


Eq. (9) takes the form 



ie „ No- E ' 


. 2 .. 2 


k G a ( k > CT ) + m k 


v a / 2 



( 11 ) 


where 


E 


Z 



ej Fj (k, o, 


s ) 


( 12 ) 


In writing Eqs. (11) and (12), we have assumed that k is along the z-axis. Fol- 
lowing the procedure outlined in references 12-14, we can easily evaluate the 
collision term and on integrating Eq. (11), we get 


F a (k, a, s) 


•/* 


ie 8 NE z ' 

«,W*Tr P.(«) 


a s cr 


kL 


: / k 


j da; e" Say k Jdn [f, (H ) (K x b + K 2 a + K 3 a ) 


+ F b 



a / b , (13) 
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where o' - (a x , a y , ct z ') , 


H. 


QAo) = T dcr 1 G (k, o * ) e 


.'A 


(14) 


P a (•) 


■-t 


/•00 

, 2/2 

r 

SCT z V a 

dcr 2 ' aj exp 

k 2 




(15) 


K b (o , n) 


27T‘ 


/ m £ 


g • O in) {g • (g ~n)}' 


b / 0~>i) (»-n) 


exp 


- ~y O - n) : 


- (16) 


K 2 a (a,n) = 


77 


g • ( g ~n) _ {g • (g ~n)} 2 
(g-n) 2 2(g-n) 4 


exp 


_ 2 


(17) 


K, a (g , n ) - 


77 


(g • n) (g • n) : 


- a 
Z7/ 


exp 


— ( g _ n ) : 


(18) 


and 


K 4 a (g, n) 




O ’ n) 


O ’ n) : 


exp 


- ~2 O _t l) : 


( 19 ) 


In Eq. (13), we have introduced the Coulomb mean free path L a - v a 4 y (Nr a ) . 
If we assume that the collisions are infrequent i.e., kL a » 1 , then to first 
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order in (kL a ) 1 , Eq. (13) gives 


F (k, a, s) = 


so- z /k 


ie a NE 2 ' 

Qa (•) + m k P a (•) 


V a Scr, /k 

klT e 


d,; Jd, e 


S7, 2 /k 


(K x b +K 2 a +K 3 a ) 


ie a E z' N I f ie b E z' N 

iQ a OO + mgk P a (n)^K 4 a jg b („) + mtk P b (n)] 


( 20 ) 


On substituting Eq. (20) in Eq. (12), we obtain 


F , _ s ) 

W(k, s) ’ 


( 21 ) 


where 


_ , . 4i TT 

4>(k, s) - + -g~ 


y Ve e f , -so- '/k r 

/ e j Qj ( °) + kL e J d<T z e " J dn e 
f 

jo. (1) (v *K‘ ^ Q, (n) (K,- +k 2 ‘ + *,■),.„ 


S7 >zA 



and 


^ — i 0) o j r v*ou 2 r 

T(k , s) = 1 - 21 IT P> <• ' o> - -W-J 

J-e, i e •'O 


- sc 

dcr e z 

z 


> 


Jdq e 


5 7J z /k 


P e (r,) (Kj 1 + K 2 e + K 3 e ) ff _ 0 - Jp, ( n )K 4 e (a = 0) 


m ^ J 

^P e (n)K 4 i (a = 0) + — Pi (»1 ) (Kj 6 + Kj + K 3 i ) ff _ ( 


(2 


In Eq. (23), ,:< p 2 j 4 tt Ne 2 /nij is the plasma frequency for the species J. We 
may note here that <£(k, s) depends only on the initial perturbation. If we con 
only those perturbations for which <h(k, S )is analytic in the complex s-plane, 
then for the Laplace inversion of Eq. (21), we have to consider only the zeros 
^(k, s) which are given by 


£ 

J=e, i 


CO 


pj 


P T (a = 0) 


3 2 

V J V 00 
e c pe 


(a . 

\ ei 


+ a. . 

1 1 


+ a 


a. ) 
ie/ 


where v c - v e jh e is the effective electron collision frequency and 


a 


■ r 

•'O 


dcr 


z'A 


fa 




P e (n)K 1 i (a = 0) 


m 

“ m p i (n)V (a - 0) , ( 
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(26) 


a ii = £ j o do-;e' Sa ’ /k Jdne S ^ /k [W+K^.JP^,) , 

a ee = | dcr z 'e" Saz ' A Jd n e S ^ /k P e (r,)(K 2 e +K3% =0 


(27) 


and 


a. 

i e 


1 00 

, -so-' /k f , S7J /k m ^ 

^ da z 'e z/ Jdne 


(o - 0) 


- P e (n)K 4 i (a = 0) • (28) 


Eq. (24) represents the required dispersion relation which gives the oscillatory 
behavior of the plasma under consideration. 


m. ANALYSIS OF THE DISPERSION RELATION 


From Eq. (15), we have 


Pj (® = 0) = 


where a - s/(kvj) and 


[1 


(29) 


Er f ( 


■>-f 

x 


dy e" 


(30) 
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While evaluating these integrals (Pj etc.), strictly speaking we should take into 
account the pole contribution which gives Landau darning but since this has been 
discussed earlier, we shall not repeat it here. We shall further assume that 
a. » 1 but a e << 1 which actually is a consequence of the fact that v e 2 » v. 2 ; 
thus on using the proper asymptotic and series expansions for the error 
function, 17 we immediately get 


Pi O = 0) 



and 


P e (» = 0) = 



(31) 


(32) 

f 


Like in the collisionless theory, we shall neglect terms of the order of a e 2 . 
Moreover, in evaluating the collision integrals, we shall retain terms only 
upto (1/a. 2 ) . 

Let us first consider the electron-ion contribution which is given by a ei in 
Eq. (25). On substituting Eqs. (16) and (19) in Eq. (25), we obtain 


a 


ei 






(33) 


10 



where 
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and on using Eq. (36) for P i5 Eq. (35) on performing ^ -integration reduces to 



r 

•'n 


dy exp 


4 {v + b 2 ) 


v 2 v 2 
sy v i y 


1 - very + J3y 


f3 f \ v cr B 3 

+ Jfr&j 1 - - 3w j - 4 (v + b 2 ) 2 


(37) 


where 


= 



and 



(38) 


Now to carry out the y-integration, we shall neglect terms of the order of 
8 2 = v i 2 /v e 2 as well as terms of the order of (l/ a . 4 ) . The resulting expression 
can then be put into the dimensionless form by using the variables (v e cr) = x 
and 1/2 v e 2 ( v + b 2 )~ 1 = z 2 ; Eq. (37) thus reads 
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where 



and 


b 


3 



(40) 


I j 1 * can be evaluated by proceeding on the lines similar to the ones for I 4 e and 

we get 



The ion-ion contribution can be similarly obtained by simplifying Eq. (26); 
this gives 


a. . 


1 l 


- 8/(l5rr 1/2 v { 6 a. s ) 


(42) 
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This, with i replaced by e, is exactly the same as the electron-electron contribu- 
tion obtained by Comisar. 12 Eq. (28) however, on integration gives 



where 



Now it is easy to see that S max is 1/S; so we immediately find that a. e ~ (a e 3 M/m) a., 
and can thus be neglected. On using the similar argument, we can show that 
a ee ^ ( a e ? “■) Nonce to ?.ri3ly7.fi the dispersion relation, we have to con- 
sider only the ion-ion and the electron-ion collisions. 

Let us now introduce the following dimensionless variables: 


CO 

co . 

p» 


- CO , 


CO . 

pi 


B, x ~ T and 


a - k/ e , 
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where ^ e 2 - KT e ^/(47rNe 2 ), is the Debye length for the electrons. Eq. (24), 
with the help of Eqs. (31), (32), (39), (41) and (42) can be written as 


1 



( 45 ) 


where 


B€ 

C = (Cii +C ei ) ’ (46) 

with e = m/M, C u = v e 6 a.. and C ei = v e 6 a ei . We shall solve Eq. (45) by the 
method of successive approximations. To lowest order, on neglecting C and a e , 
we obtain 


c °* 2 = [i + (1. + 12 v a2 T ) 1/2 ] s % 2 say • (47) 

Here y stands for (l + 1/a 2 ). For T « 1, Eq. (47) further reduces to 

w o 2 = i + 3k2 ^i 2 : (48) 

the result obtained by Bernstein and Trehan. 18 k i in Eq. (48) is the ion-Debye 
length. To next order Eq. (45) takes the form 
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which we shall solve numerically. It is worth pointing out that 


Cu(“* = "o) = - 


8i a 5 

15tt 1/2 T 1/2 ^ 5 


( 50 ) 


being pure imaginary, does not affect the real part of co ; however, C ei = ^ 0 ) 
is a complex quantity and thus changes both the real and the imaginary parts of 
co*. C.j and C ei for various values of a and T are given in Table 1. 

According to Eq. (49), although Re co [co t ) increases with an increase in a, 

T as well as B but this increase, as shown in Tables (2a) and (2b), is negligible. 

As illustrated in Figures 1 and 2, Im •■co* {co*) increases with increase in a and 
B but decreases with increase in T. So two-body Coulomb collisions have a tendency 
to stabilize the ion acoustic waves. As T increases, Landau damping takes over 
the collisional damping. For the sake of completeness, in Figure 3, we have 
shown the thermal effects on the ion waves. 

In a system where the collisions are frequent, this model will break down; 
in this case one should use the Kinetic equation which tak^s into account the 
correlations between the charged particles. 


IV. CONCLUSIONS 

Independent of the ratio of the ion temperature to the electron tempera- 
ture T, the characteristic frequency of the ion acoustic waves in a plasma with 
weak Coulomb collisions gets affected only by the electron-ion collisions. 
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Just as in a collisionless plasma, in this case also a> T increases with increase 
in T; moreover it increases with increase in the collision frequency v c as well 
as with (kX. e ) . The electron-electron and the ion-electron collisions are 
negligible compared to the ion-ion and the electron-ion collisions. The latter 
two play a somewhat equally important role in damping these waves. The col- 
lisional damping increases with increase in v c and (k\ e ) but decreases with 
increase in T. 
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Table 1 


C u andC ei for Various Values of a - kk e and T - T i/ /T e 
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Table 2a 


Rea;* for a = 0.1 and T = 0.1, 0.05 and 0.01 for Various Values of B 


B 

T = 0. 1 

T = 0.05 

T = 0.01 

0.0 

0.110936 

0.105914 

0.100927 

0.0001 

0.110938 

0.105916 

0.100929 

0.001 

0.110959 

0.105939 

0.100953 

0.01 

0.111167 

0.106160 

0.101187 

0.05 

0.112108 

0.107161 

0.102213 

0.10 

0.113323 

0.108447 

0.103456 


Table 2b 


Re co for 

a = 1.0 and T = 0.1, 0.05 and 0.01 for Various Values of B 

B 

T = 0. 1 

T = 0.05 

T = 0.1 

0.0 

0.883308 

0.858313 

0.838316 

0.0001 

0.883309 

0.858313 

0.838317 

0.001 

0.883315 

0.858318 

0.838322 

0.01 

0.883376 

0.858370 

0.838368 

0.05 

0.883647 

0.858599 

0.838566 

0.10 

0.883987 

0.858881 

0.838807 
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ILLUSTRATIONS 


Figure 1. 


Figure 2. 


Figure 3. 


Variation of Im a? with B for T = 0.1 is illustrated for a = 0.01, 0.1 
and 1.0 

Variation of Imw withB for a = o.l is illustrated for T = 0.01, 0.05 
and 0.1 

Variation of Re w with a for B = 0 is illustrated for T = 0.01 and 0.1 
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